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Abstract 

We study a continuous-space version of the totally asymmetric simple exclusion process 
(TASEP), consisting of interacting Brownian particles subject to a driving force in a 
periodic external potential. Particles are inserted at the leftmost site at rate a, hop to 
the right at unit rate, and are removed at the rightmost site at rate /3. Our study is 
motivated by recent experiments on colloidal particles in optical tweezer arrays. The 
external potential is of the form generated by such an array. Particles spend most of 
the time near potential minima, approximating the situation in the lattice gas; a short- 
range repulsive interaction prevents two particles from occupying the same potential well. 
A constant driving force, representing Stokes drag on particles suspended in a moving 
fluid, leads to biased motion. Our results for the density profile and current, obtained 
via numerical integration of the Langevin equation and dynamic Monte Carlo simulations, 
indicate that the continuous-space model exhibits phase transitions analogous to those 
observed in the lattice model. The correspondence is not exact, however, due to the lack 
of particle-hole symmetry in our model. 
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I. INTRODUCTION 



A driven lattice gas, or driven diffusive system, is a system of interacting particles 
that jump in a preferred direction on a lattice. The system cannot reach equilibrium 
but may attain a stationary state with a steady current; the model is a prototype 
for studies of nonequilibrium states [l-4]. The simplest example of a driven diffusive 
system, which has become one of the standard models of nonequilibrium statistical 
mechanics, is the totally asymmetric simple exclusion process (TASEP) 5|-[8(. In the 
TASEP with open boundaries the edge sites are connected toparticle reservoirs with 
fixed densities. Introduced as a model of biopolymerization [9| and transport across 



membranes JlO| , over the years, this mode 
traffic flow and cellular transport 
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ras been applied to other processes, e.g., 



From a mathematical point of view the model is of interest in the theory of in- 
teracting particle systems since, despite its simplicity, it shows a nontrivial behavior 
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15| . In the one-dimensional TASEP with open boundaries, particles jump 



only to the right, along a one-dimensional lattice whose sites can be empty or occu- 
pied by a single particle. Particles are injected at the leftmost site at rate a if this 
site is empty, and removed at the rightmost site at rate /3 if this site is occupied. 
The one-dimensional TASEP, which has been solved exactly, exhibits three distinct 
phases in the plane 
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19|. The phase transition is discontinuous along the 



line a = (3 < 1/2, where the density profile is linear, and continuous along the 
lines a = 1/2, a > (5 and /3 = 1/2, a < (3 (see Fig. |9]). Although this model and 
and variants have been the subject of intensive theoretical study, there is as yet no 
realization of a TASEP-like system in the laboratory. 

The invention of optical tweezer arrays has permitted investi gati on of the dy- 



namics of colloidal particles in an external periodic potential 



20 



-|25|. The motion 



at long times and low friction consists of jumps between adjacent potential min- 
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ima. If the particle and potential well sizes are chosen properly, only one particle 
can occupy a given well. The exclusion process is a caricature of this dynamics, 
suggesting that a system of colloidal particles in an optical tweezer array could be 
designed as a laboratory realization of the TASEP. Motivated by this possibility, 
we propose a model in continuous space having the same essential characteristics 
as the lattice TASEP. Our model represents colloidal particles immersed in a fluid 
flowing at constant rate through a one- dimensional optical tweezer array, restricting 
particle motion to the array axis. Specifically, we study a one-dimensional system 
of interacting Brownian particles subject to a periodic external potential, and to a 
constant external force representing the drag due to the fluid motion. A short-range 
(essentially hard-sphere) repulsion between particles prevents more than one parti- 
cle occupying the same potential well. This continuous-space model is studied via 
numerical integration of Langevin equation and dynamic Monte Carlo simulation. 
We observe phase transitions similar to those found in the lattice TASEP. Some 
differences in the detailed behavior nevertheless appear, due to the lack of particle- 
hole symmetry in the continuous-space model. Details on the model and simulation 
methods are given in Sec. [Til and IHH Simulation results are presented in Sec. IIV| 
while our conclusions and prospects for future work are outlined in Sec. |V] 



II. CONTINUOUS-SPACE MODEL 



Our aim is to study a continuous-space model sharing the same essential features 
as the TASEP (defined on a lattice), as a first step toward experimental realization 
of a TASEP-like system. The model should possess the following characteristics: i) 
confinement of particles to a one-dimensional structure; ii) localization of particles 
at potential minima ( "wells" ) of a linear periodic array with iii) multiple occupancy 
prohibited; iv) biased hopping between adjacent wells; v) insertion (removal) of par- 
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tides at the initial (final) well. Criteria i)-iii) are realizable with a suitably tailored 
optical tweezer array. The array consists of a series of spherically symmetric op- 
tical tweezers; the particles flow along this line, which we take as the x axis. To 
avoid particles escaping the array, there should be a substantial overlap between 
neighboring wells, so that a potential maximum at a point midway between two 



wells is in fact a saddle point in the full three-dimensional space (see Fig. 1(a)). 
For a TASEP-like system, it is crucial that the probability of a particle escape from 
the array be negligible on the time-scale of the experiment. Effective confinement 
requires a potential barrier to escape the array much larger than ksT; the barrier 
between adjacent minima should be smaller, to allow transitions between neighbor- 
ing wells. In what follows we shall assume this condition is satisfied, and consider, 
for simplicity, a one-dimensional system. Fluctuations of the particle positions in 
the directions perpendicular to the array will therefore be ignored, but should be 
included in a more complete analysis. 

The diameter of the optical tweezer well should be slightly greater than the par- 
ticle diameter, so that at most one particle can occupy the well at a given time. 
Due to thermal fluctuations, particles can occasionally overcome the potential bar- 
rier separating neighboring wells. Since hopping must be asymmetric, we impose 
a steady fluid motion along the +x axis, which effectively prohibits particle jumps 
in the opposite direction. In the lattice TASEP particles are inserted in the first 
site and removed from the last. Experimental realization of this feature is subtler, 
but can in principle be achieved with the help of optical tweezers at the beginning 
and end of the array, which drag particles into the first well and out of the last one 
at prescribed rates. We discuss an alternative method of insertion and removal in 
Sec. El 

The above sketch of an experimental setup motivates our study of Brownian 
motion of interacting colloidal particles in an optical tweezer array. The Langevin 
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equation for the i-th particle is 



mxi(t) = -b [±i(t) -v\- dVe ^ x ^ _ JL [V int (x iti -i) + V int (x i4+ i)] + mTi(t), (1) 

OXi OXi 

where Xi, Xi and ii are, respectively, the position, velocity and acceleration of par- 
ticle z, and v is the terminal velocity of a particle in the moving fluid, 
in the absence of the periodic external potential V ext . Vi nt is the (strongly repulsive) 
potential between neighboring particles. (We assume that the range of Vi n t to be 
short enough that only neighboring particles interact.) The first term on the right 
side of Eq. ([T]) represents damping of the particle velocity relative to the fluid. For 
a sphere of radius R, Stokes' Law gives: 

b = 6tt/iR, (2) 

where /i is the fluid viscosity. Here it is important to stress that the fluid is three 
dimensional although we treat the particle motion as one dimensional. The final 
term is a random noise with the following properties: 

(TO = 0, (3) 

(r^r^O) = ^^M(*-0. ( 4 ) 

where ks = 1.3806504 x 10~ 23 J/K is Boltzmann's constant, T is temperature, and 
m is the particle mass. It is convenient to ignore the inertial term mxi(t) in Eq. (pQ), 
since the observational times of interest (microseconds or greater) are much larger 
than the relaxation time of the velocity, m/b ~ 10~ 8 s, for our choice of parameters. 
Then the velocity of particle % follows, 

Xi{t) = v - ^ ~ \~^T l V int{xi,i-i) + V ir a(x it i + i)] + y r <(*)- ( 5 ) 
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A. Potentials 



The potential of an optical tweezer array can be represented by a sum of n 
identical Gaussian profiles of width a and spatial period d: 

N 

-(x-dn) 2 /2a 2 

n=0 



V t (x) = -V J2e- 



(6) 



We require neighboring wells to overlap, which can be accomplished by setting a ~ d; 



Figure 1(a) shows that this condition is satisfied for a = d/A. For this choice of 
parameters, the potential is well approximated by a cosine, as can be seen from the 
Fourier coefficients 

a 



a, 



\ 1 



N 
n=0 



-{x-dn) 2 /2a 2 



dx. 



(7) 




FIG. 1: (a) Two neighboring Gaussian profiles (red and black curves) and their sum 
(green); (b) effective external potential. 

Numerical evaluation yields 
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a = -0.627 



ai = 0.365 



(8) 



a 2 = 9.010 x 10~ 3 . 



Note that 02/0-1 — 0.025, allowing us to write, to a good approximation, the array 
potential as a constant plus a cosine term. For the one-dimensional model studied 
here, we define 



as the periodic external potential. 

The interaction between colloidal particles is taken as purely repulsive; for con- 
venience we use a truncated 1/r 12 potential: 



where a is the particle diameter and r distance between neighboring particles. 
B. Parameter Values 

To specify the external and interaction potentials, we need to fix Vo and Uq. 
These values must be chosen so as to approximate the TASEP dynamics given the 
length, time and energy scales characterizing the system. We measure lengths in 
units of microns, time in seconds and energies in units of k^T = 4.1419512 x 10 _21 J, 
assuming a temperature of 300K. We set d = 1 (so that the period of the external 
potential is 2 microns), and take the particle diameter as a = 1.8, so that a pair 
of particles occupying neighboring wells have some freedom to fluctuate about the 
potential minimum. Taking the fluid as water (with viscosity \x = 0.01g/cm-s) we 
fix the friction coefficient as b = 4.096 A; j gT-s/(yum) 2 . 




(9) 




(10) 
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To determine the fluid velocity v and external-potential intensity Vq we examine 
the effective external potential, defined as 

i.e., the sum of external periodic potential and a fictitious potential representing the 



constant friction force acting on particles (see Fig. 1(b)). If bvdo S> Vb, particles 
do not feel the periodic potential, while if bvdo <C Vq particle hopping is essentially 
unbiased. Let AEq (AEp) denote the difference between a given maximum of V e ff 
and the first minimum to the right (left) of this maximum. In this way, AEq (AEp) 
is the potential barrier separating a given potential minimum from its left (right) 
neighbor. A simple calculation shows that 

AE G - AEp = 2bvd , 

(12) 

AE G + AEp = AVo cos(7rx /rfo) + 46ra , 



d . _! / bvdo , 
x = — sin — — , (13) 



with 

Xo = 

TX \ TTVo 

such that the potential minima occur at Xj = 2jdo + xo, and maxima at yj = 
2jdo + 1 — xq for j an integer. Given AE G and AEp, we can determine xq, v 
and Vq. We take AEp ~ kpT, and AE G ^> kpT, so that a particle has a finite 
rate of jumping to the well on the right, and virtually no chance of jumping in the 
opposite direction. A good correspondence with lattice TASEP is obtained using 
AEq = 20kpT and AEp = 2fc#T. For these values we have, 

v = 2.1974 jum/s 

V = 4.5680 kpT (14) 
xo = 0.2158 fim. 

To maintain the interparticle repulsion in the presence of the external potential, we 
must take Uo substantially greater than Vq. On the other hand, very large values 
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of Uq are inconvenient for numerical integration of the Langevin equation, as a very 
small time increment would be required to avoid spurious particle displacements. 
We therefore use U = 40k B T ~ 9V . 

C. Single-particle dynamics 

To begin, we consider a single Brownian particle moving in the system. Figure [2] 
shows a typical evolution of the particle position over time. The graph exhibits 
plateaux whose size corresponds to the time a particle stays in a given well. Studies 
of this kind allow us to determine the mean transition time r between neighboring 
wells as 6.5380(9)s. This quantity is needed in order to define the insertion and 
removal rates. (Recall that in the lattice TASEP these rates are defined in units 
of the hopping rate). If the first well is empty we insert a particle there (at the 
potential minimum position), at rate a/r; if the last well is occupied, we remove 
the particle at rat 




Time [s] 



FIG. 2: Typical particle trajectory. 
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III. DYNAMIC MONTE CARLO SIMULATIONS 



The Langevin simulation outlined in the preceding section is valuable for fixing 
the time scale of hopping between wells, and for confirming the basic phenomenology 
of the model. It is, however, rather inefficient numerically, so that it is desirable 
to implement a dynamic Monte Carlo (MC) simulation for large-scale studies. We 
apply the Metropolis algorithm to evolve the particle positions in time. 

In this approach we use the following expression for the potential energy: 

N N-1 
E = -^[Vo COs(7TXi) + bvXi] + ^ V int (Xi + i - Xj), (15) 

1=1 1=1 
where V, n t is given by Eq. ([TO]) . The drag force due to the moving fluid is repre- 
sented by the effective potential —bvx. In the MC dynamics, a trial configuration 
is generated by selecting one of the N particles at random and subjecting it to a 
random displacement Ax, chosen from a Gaussian distribution with mean zero and 
standard deviation a = 0.2/im. This value is 10% of the well size, large enough to 
afford a substantial speedup, but small enough that the probability of a particle dis- 
placement greater than 2/zm is negligible. As is usual in Metropolis MC, trial moves 
such that the change in energy AE < are always accepted, while for AE > the 
trial move is accepted with probability e - A£; / fe s T . We determine the mean number 
of Monte Carlo steps required for a particle move from one well to its neighbor on 
the right as N MC = 120.917(1). Thus the time per Monte Carlo step is 

T = — — = 5.4070(8) x 10~ 2 s. (16) 

Nmc 

Particle insertion and removal are done as in the Langevin simulations. We verify 
below that this method is equivalent to the latter approach; the MC algorithm is 
~ 1000 faster than numerical integration of the Langevin equation, for L=100. 
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FIG. 3: Density profile for (3 = 0.4 obtained via numerical integration of the Langevin 
equation, for a values as indicated, L = 100. 

IV. RESULTS 

During the simulations we monitor the mean occupation probability p(i) at each 
well i, and the current J, measured by the mean number of particles leaving the 
system per unit time. Examples of density profiles in the stationary regime {p{i) 
versus i), are shown in Fig. [3] for several values of a, with /3 = 0.4. (From here on 
a and /3 are given in units of 1/r, where r is the mean time required for hopping 
between wells.) This figure shows that the continuous-space model exhibits the same 
basic phenomenology as the lattice TASEP. For a < 0.4 the overall density grows 
with a. On increasing a from 0.4 to 0.5 there is a marked increase in density, but 
for further increases the density changes very little. While the Langevin simulation 
(LS) results already suggest that the model exhibits phase transitions, we shall use 
the more precise results of our MC simulations to perform a detailed analysis. Before 
proceeding, we verify that the MC method yields results in agreement with the LS. 
In Fig. H] we compare density profiles obtained via LS and MC for the same values 
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of a and (3. The bulk densities obtained using the two methods differ by < 1.6%. 
Thus the MC method captures the behavior found using the Langevin equation to 
good precision. 
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FIG. 4: Comparison between profiles obtained using LS (open symbols) and MC (filled 
symbols), a = 0.3 (lower curves), a = 0.8 (upper curves), /3 = 0.4 , L = 100. 

We perform MC simulations of systems of L = 100, 200 and 500 wells. Far from 
the phase transition, density profiles depend only weakly on system size, but near 
the transition there are significant finite-size effects, as illustrated in Fig. [5j In this 
case, as L increases, the profile tends to a near-constant value except for a sharp 
increase near the exit. 



To minimize boundary effects we study the bulk density p, defined as the mean 
density over the 10% of sites nearest the center: 

1 „ 0.55L 

p = { E *>(«')• 

i=0.46L 
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FIG. 5: Density j id 500 wells. 
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a 

FIG. 6: Bulk density versus a for (3 = 0.4, system sizes L = 100, 200 and 500. 

The bulk density as a function of a, for = 0.4, is shown in Fig. [6] for the three sys- 
tem sizes studied. These results strongly suggest the development of a discontinuity 
in p(a) near a = 0.47 as the system size is increased. 
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(c) (d) 

FIG. 7: Upper panels: bulk density for (a) fixed a and (b) fixed /?. Lower panels: current 
for (c) fixed a and (d) fixed (3. System size L = 200. 

A. Phase diagram 

Of principal interest in determining the phase diagram are the bulk density and 
the current as functions of the rates a and /3. These results are summarized in 
Fig. [7J showing evidence of both continuous and discontinuous phase transitions, 
depending on the rates. We see that for low a (a < 0.8 or so) and (3 < a the system 
is in the high-density phase, in which density and current depend only on /3, whereas 
for /3 > a, and j3 < 0.6 or so, the system is in the low-density phase in which density 
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and current depend only on a. For larger values (a > 0.8 and (3 > 0.6), the system 
is in the maximum-current phase, in which density and current are independent of 
both a and (3, and the current takes its maximum value. Thus the continuous-space 
model exhibits the same three phases observed in lattice model. As in the lattice 
model, the transition between the low- and high-density phases is discontinuous, 
whereas transitions between the maximum-current phase and the other phases are 
continuous. (While the density is discontinuous in the former case, the current is 
always continuous at the transition.) 

We adopted the following procedure ("polynomial method") to determine the 
values of (a, /3) along the discontinuous transition line. Consider the case of fixed a. 



For small /3, the stationary current depends only on (3 (see Fig. 7(c)). We therefore 
fit a polynomial P{(3) to the current in this regime, using data for a large, fixed a 
(in practice, a > 1.4). For larger values of j3, the current depends only on a; in 



Fig. 7(c) this regime corresponds to one of the plateaux, J = J* (a). The transition 
point /3 c (a) is taken as the value at which the plateau intersects the polynomial fit to 
the small-/? data, i.e., P((3 C ) = J* (a). Determination of a c (f3) follows an analogous 
procedure, in which we fit a polynomial to the current data for small a. We find 
that a quadratic polynomial is sufficient to fit (to within uncertainty) the current 
j(j3) at fixed a, while a good fit of j(a) (at fixed 0) requires a quartic polynomial. 
In both cases the constant term in the polynomial is zero, since the current vanishes 
for a and/or /3 zero. 

The above method is quite effective in locating points along the discontinuous 
transition line. Although it can in principle be used to locate continuous transitions 
as well, we found that in this case the estimates for a c and /3 C are rather sensitive to 
one's choice of the range of values fit using the polynomial. We found the following 
approach ("derivative method") to be more useful for continuous transition points. 
For fixed a, we estimate the derivative dJ/d(3 using a spline fit (see Fig. [8]). The 



15 



derivative decreases in a linear fashion with increasing /3, except for a small roundoff 
region that we interpret as a finite-size effect. We then estimate (3 C as the point where 
dJ/dP falls to zero, using the data in the linear region. The procedure for fixed j3 
is analogous. The continuous transition points obtained by this procedure are quite 
robust with respect to changes in the region analyzed, as long as we exclude the 
roundoff region. 

0.15 
0.1 

ca 

— ! 

0.05 



w () 0.5 1 1.5 

P 

FIG. 8: Current derivative for a = 2.0 and L = 200. 

Using the methods described above we construct a phase diagram based on the 
data for each system size studied; that for L = 500 is shown in Fig. The vertical 
boundary at a = 0.905(2), and the horizontal boundary at /3 = 0.714(1) are obtained 
using the derivative method. (The polynomial method yields a = 0.8(1) and j3 = 
0.67(5), respectively, for these boundaries.) The grey triangle in Fig. |9] represents a 
region on which we could not determine precisely the phase boundaries using either 
method, due to numerical uncertainty and finite-size effects. While the simplest 
interpretation is that the three phase boundaries meet at the point of intersection 
between the continuous transitions (i.e., horizontal and vertical) lines, the data in 
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FIG. 9: Phase diagram of the TASEP on the lattice (dotted lines) and in continuous- 
space (solid lines) for L=500. The grey triangle denotes a region in which we are unable 
to determine the phase boundaries precisely. 

hand are not sufficient to verify this. 

The phase boundaries for L = 100 and L = 200 are nearly the same as for 
L = 500, suggesting that the latter are already quite close to limiting (infinite-!/) 
values. Table [J gives the a and (3 values for the continuous transition lines for the 
three sizes. There is little sign of systematic variation with system size, other than 
a small (~ 2.5%) increase in j3 c on going from L = 200 to L — 500. 

While the continuous-space phase diagram is isomorphic to that of the lattice 
model, there are some differences between the two cases that appear likely to persist 
in the infinite-size limit. Since the lattice model possesses particle- hole symmetry, 
the phase diagram is invariant under the exchange of a and /3. Thus the boundary 
between the high- and low-density phases is a straight line extending from the origin 
to the point (1/2, 1/2) in the a-/3 plane. The phase diagram of the continuous-space 
model does not possess this symmetry; the phase boundary between the high- and 
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100 0.901(4) 0.701(2) 
200 0.901(2) 0.697(1) 
500 0.905(2) 0.714(1) 

TABLE I: Values of a and ft at the continuous transitions obtained via the derivative 
method. 

low-density phases does not fall along the line a — ft, and appears to be somewhat 
curved. 

One might inquire whether the differences in the phase boundaries of the lattice 
and continuous-space models merely reflect finite-size effects in the latter. We have 
verified that in the lattice TASEP with L = 200, the phase boundaries fall quite 
near their expected (infinite-size) positions. Comparison of the phase boundaries (in 
continuous space) for L = 200 and 500 suggests that finite-size effects are somewhat 
stronger in the continuous space model than on the lattice. Given the lack of particle- 
hole symmetry, however, it appears very unlikely that the continuous-space phase 
boundaries will converge to those of the lattice model in the infinite-size limit. 

The differences between the lattice and continuous-space models reflect, in part, 
the absence of particle-hole symmetry in the latter; particle positions fluctuate in 
continuous space, but are fixed in the lattice model. In continuous-space, moreover, 
particles occupying neighboring wells may influence one another via the repulsive 
potential Vint- On the lattice model no such influence exists, beyond simple exclu- 
sion. In continuous space, repulsive interactions should tend to spread particles more 
uniformly than on the lattice, promoting particle removal, and hindering insertion. 
Thus the transition from high to low density occurs for ft < a. Since repulsion is 
more significant for higher densities (i.e., larger a) the phase boundary should curve 
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toward the a axis, as is observed. The smaller value of (3 at the continuous transition 
(approximately 0.714(1) for L = 500), compared to that of a (about 0.905(2) for 
L = 500) may also be attributed to repulsion between neighboring particles. 



B. Current: comparison with mean-field theory 

The absence of particle-hole symmetry is again evident in a plot of the current 
as a function of density. In the lattice model, mean field theory gives J = p(l — p) 
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271 ] which is in fact an exact expression. Fig. [10] compares the current on lattice 
with our results for continuous space. (Note that the latter exhibit virtually no 
finite-size effects on the scale of the figure.) Unlike in the lattice TASEP, here the 
current is not symmetric about p = 1/2; it takes its maximum value at a density of 
about 0.57. The fact that the maximum current occurs at a higher density than on 
the lattice may again be attributed to interparticle repulsion. 




FIG. 10: Current versus density in the lattice model (exact) and the continuous-space 
model. 
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On the lattice, the current is equal to the probability of having an occupied site 
with its neighbor on the right vacant, i. e., J = (£t(l — £t+i))> Vi, where & is an 
indicator variable equal to one if site i is occupied and if it is empty. In mean-field 
theory the joint probability is factored so: — = — £t+i)), and 

setting = p in the bulk, we obtain J = p(l — p). 

In developing a mean-field theory for the continuous-space model, one might 
argue that J ~ p(l — p)/r, i. e., that the current is simply given by the transition 
rate 1/r for jumps between neighboring wells times the probability of a given well 
being occupied and its right neighbor empty. Due to repulsive interactions between 
neighboring particles, however, a particle in the well to the left can also influence 
the hopping rate. As a first approximation we write 

J = P(0, 1, 0)j(0, 1, 0) + P(l, 1, 0)j(l, 1, 0), (18) 

where P(£j_i, is the joint probability for three adjacent wells, and 
is the transition rate in this configuration. Factorizing the joint 
probability, we have P(0, 1, 0) = p(l - p) 2 and P(l, 1, 0) = p 2 (l - p). It remains to 
evaluate the currents j(0, 1, 0) and 1, 0). 

Consider first j(0, 1,0), the rate to overcome the barrier between wells % and i + 
given that both i — 1 and i + 1 are empty. The mean first-passage time r ca , for a 
particle to overcome the barrier, is readily found via analysis of the one-dimensional 



Fokker-Planck equation. From the standard result 



28j we have 



U{ X >)/k B T dx , / e -U{x")/k B T dx ,^ (lg) 



k B T 

where a = Xq and c = Xq + 2d are the positions of adjacent potential minima, 
with Xq given by Eq. f TTS]) . Using the effective external potential, Eq. (II ip . for 
U(x) (since there are no interactions with other particles), numerical evaluation of 
Eq. (1191) yields r ca = 6.538 s= l/j(0, 1,0). This value agrees to within uncertainty 
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with our simulation result for the mean time r for a particle to hop between adjacent 
wells when there are no other particles in the system. Thus the mean-field curve in 
Fig. [10] agrees with simulation in the low-density limit. 

To estimate the transition rate j(l, 1,0), we perform a Monte Carlo simulation 
to determine the mean time required for a particle to hop to the next well, when 
the preceding well is occupied, and there are no other particles in the system. This 
yields l/j(l, 1,0) = 5.548s. Although the presence of the trailing particle leads 
to an increase of about 18% in j, the effect is not sufficient to yield quantitative 
agreement with the current observed at higher densities (Fig. [TO]) . There are two 
possible sources for this discrepancy. First, at moderate and high densities, strings 
of n > 3 occupied wells occur with finite probability, and the cumulative effect of 
repulsions along the chain should make the hopping rate of the lead particle an 
increasing function of n. Simulations of n = 3 - 9 occupied wells show that the 
transition rate of the first particle grows with n, but not enough to account for the 
maximum value of the current observed. 

A second point is that the mean-field factorizations P(0, 1,0) = p(l — p) 2 and 
P(l, 1,0) = p 2 (l — p), are not very accurate in the continuous-space model. For 
density p ~ 0.52, for example, we find 

P(0, 1,0) . , 

1.21 (20) 



p(l-p) 2 
and 

P(0,1,1) 



1.04 (21) 



p 2 (l-p) 

implying a significant correction to the mean-field theory predictions. 
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C. Fluctuations 

To close this section we note an interesting finding on fluctuations. The variance 
of the density, as a function of a, with fixed /3, exhibits a maximum at the discon- 
tinuous transition (see Fig. [TT]) . For values such that the transition is continuous, 
by contrast, no peak in var(p) is observed at the transition. (Similar behavior is 
found, varying with a fixed.) The large density fluctuations are associated with 
the presence of a shock separating high- and low-density regions, whose position 
fluctuates over the entire system. The position of maximum variance agrees to 
within uncertainty with the lines of the (discontinuous) phase transitions reported 
in Fig. [9j Analysis of the variance, however, appears to furnish less precise results 
than the method described above. The total energy exhibits fluctuations similar to 
those observed in the density, but we do not find any signal in var(J) associated 
with the phase transitions. 
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FIG. 11: Density variance versus a for (3 values as indicated. System size L = 200. 
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V. DISCUSSION 



We propose a continuous-space model of interacting Brownian particles in a pe- 
riodic potential as a possible realization of the TASEP. The particles are subject 
to a constant drive in a periodic external potential. Using numerical integration 
of the Langevin equation and Monte Carlo simulation, we study systems of L = 
100, 200, and 500 wells. Our results show that the continuous-space model exhibits 
continuous and discontinuous phase transitions analogous to those observed in the 
lattice TASEP. The phase diagram of the continuous-space model is similar to that 
of the lattice model, but exhibits some differences due to the absence of particle- 
hole symmetry. This difference appears to be associated with fluctuations of particle 
positions around potential minima. Such fluctuations, together with the repulsive 
interactions between neighboring particles, cause the current to attain its maximum 
value at a density somewhat greater than 1/2, the density marking the maximum 
current in the lattice model. We expect these changes (relative to the lattice model) 
to be generic for continuous-space systems exhibiting TASEP-like phase transitions. 

We believe that the present study demonstrates the possibility of observing 
TASEP-like behavior in laboratory experiments on systems of interacting colloidal 
particles in a one-dimensional optical tweezer array. The essential features of the 
TASEP - localization of particles in potential wells, with multiple occupancy prohib- 
ited, and bias hopping along the line - are readily accomplished with a appropriate 
choice of particle, fluid, and tweezer array parameters. It is however less obvious 
how to implement random insertion and removal of particles at the first and last 
wells of the array. Particle manipulation can be accomplished using optical tweezers 
to transfer particles between wells and reservoirs. To transfer particles in a random 
fashion, these tweezers would have to be intrinsically noisy or chaotic, controlled by 
a random number generator, or driven by a noise signal. A simpler alternative may 
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be periodic insertion and removal. In this case, one inserts a particle into the first 
well (when empty) at intervals of t/oc, and checks for occupancy of the final well at 
intervals of t//3, removing the particle if the well is occupied. 

A preliminary study of the continuous-space model using periodic insertion and 
removal confirms that the three TASEP phases are again found. The density profiles 
under periodic and random particle transfer are very similar in the maximum current 
phase, where the density profile is insensitive to small changes in the insertion and 
removal rates. Small systematic differences do however appear in the other phases, 
as shown in Fig. [121 

0.96 



0.93 
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- random 
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high density phase 
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FIG. 12: Comparison of density profiles obtained using periodic (black) and random (red) 
insertion and removal, L = 100. 

Although our study strongly suggests the feasibility of a laboratory realization of 
the TASEP, a number of additional features would have to be included in the model, 
before a quantitative comparison with experiment could be made. The principal 
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modifications we expect to be necessary are study of a three-dimensional model, 
allowing fluctuations in directions perpendicular to the array axis, and inclusion of 
hydrodynamic interactions between the particles. We defer these tasks to future 
work. 
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